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Abstract 

Using data accumulated with the CLEO detector corresponding to an integrated luminosity of 
£=5.46 pb~^ on the peak of the ijj{2S) and 20.46 pb~^ at ^/s=3.67 GeV, we report preliminary 
branching fraction measurements for seven new decay modes of the tp{2S) (jjStt, t/'Stt, 2(K'^K~), 
ppK~^K^ , AAvr+vr^, ApK~^, and Api^+vr+vr") and more precise measurements of nine previ- 
ously measured modes (2(7r"'"7r~), 2(7r"''7r~)7r'', lutt^tt^, K^tt^tt^ , (pn^TT^, ujK~^K~, cpK^K^, 
ppir^TT^, and AA). We also include a study of upp and obtain an improved upper limit for cppp. 
Results are compared, where possible, with the corresponding J/ijj branching ratios to test the 
12% rule. 



*Submitted to the International Conference on High Energy Physics, August 2004, Beijing 
^On leave of absence from University of Chicago. 
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In perturbative QCD the states J /if) and il>{'2S) are non-relativistic bound states of a 
charm and an anti-charm quark. The decays of these states are expected to be dominated 
by the annihilation of the constituent cc into three gluons. The partial width for the decays 
into an exclusive hadronic state, /i, is expected to be proportional to the square of the cc 
wave function overlap at the origin, which is well determined from the leptonic width P|. 
Since the strong coupling constant, a^, is not very different at the J/tp and ip{2S) masses, 
it is expected that for any state h the J /if) and ip{2S) branching ratios are related by 



where B denotes a branching fraction, /i is a particular hadronic final state, and the leptonic 
branching fractions are taken from the PDG [1]. This relation is sometimes called "the 12% 
rule". Modest deviations from the rule are expected Although the rule works well for 
some specific decay modes of the 'ip{2S), it fails spectacularly for ip{2S) decays to final states 
consisting of one vector and one pseudoscalar meson, such as pvr. 

Values of Qh have been measured for a wide variety of final states [H 0, 0]- Most re- 
cently CLEO has measured additional 1~0~ final states at the ipi^S) and the continuum 
including the first observation of ipi^S) pn. A recent review j3| of relevant theory and 
experiment concludes that current theoretical explanations are unsatisfactory. Clearly more 
experimental results are desirable. This paper presents measurements of the following new 
decay modes of the il){2S): tjStt, t^'Stt, 2{K~^K~), ppK^K~, AAvr+vr", ApK'^, ApK^7c~^7r~ , 
and more precise measurements of these previously measured modes: 2(71+71"), 2(7r+7r^)7r°, 
uJTi^Ti^, K^K^TT^TT^ , (f)7i^7i^, ujK^K~, (f)K~^K~, ppTr^iT~ , and AA. We also measure utpp 
and obtain an improved upper limit for (ppp. Where applicable, the inclusion of charge 
conjugate modes is implied. Eleven of the modes we study have been previously observed 
at the J/ip. 

The data sample used in this analysis is obtained at the ip{2S) and the nearby continuum 
in e"'"e~ collisions produced by the Cornell Electron Storage Ring (CESR) and acquired 
with the CLEO detector. The CLEO III detector features a solid angle coverage for 
charged and neutral particles of 93%. The charged particle tracking system, operating in 
a 1.0 T magnetic field along the beam axis, achieves a momentum resolution of ~0.6% at 
p = 1 GeV/c. The calorimeter attains a photon energy resolution of 2.2% at = 1 GeV and 
5% at 100 MeV. Two particle identification systems, one based on energy loss {dE/dx) in the 
drift chamber and the other a ring-imaging Cherenkov (RICH) detector, are used together 
to separate from ir^. The combined dE/dx-KlCH particle identification procedure has 
efficiencies exceeding 90% and misidentification rates below 5% for both ir^ and K^. 

Half of the ip{2S) data and all the continuum data were taken after a transition to the 
CLEO-c detector configuration, in which CLEO Hi's silicon-strip vertex detector was 
replaced with a six-layer all-stereo drift chamber. The two detector configurations also 
correspond to different accelerator lattices: the former with a single wiggler magnet and a 
center-of-mass energy spread of 1.5 MeV, the latter (CESR-c jsf) with six wiggler magnets 
and an energy spread of 2.3 MeV. 

The integrated luminosity (C) of the datasets was measured using 77 events in the 
calorimeter [9|. Event counts were normalized with a Monte Carlo (MC) simulation based 
on the Babayaga ^3 event generator combined with GEANT-based jT]| detector modeling. 
The datasets consist of £=5.46 pb^^ on the peak of the ip{2S) (2.57 pb~^ for CLEO HI, 
2.89 pb-i for CLEO-c) and 20.46 pb^^ at v^=3.67 GeV (all CLEO-c). The nominal scale 
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factor used to normalize continuum yields to ^ji^S) data is /nom = 0.2645 ±0.004, and is de- 
termined from the integrated luminosities of the data sets corrected for the 1/s dependence 
of the cross section, where the error is from the relative luminosity uncertainty. The actual 
/ used for each mode also corrects for the small differences in efficiency between the ip{2S) 
and continuum data samples. 

Standard requirements are used to select charged particles reconstructed in the tracking 
system and photon candidates in the Csl calorimeter. We require tracks of charged particles 
to have momenta p > 100 MeV and to satisfy | cos^| < 0.90, where 6 is the polar angle with 
respect to the e"*" direction. Each photon candidate satisfies > 30 MeV and is more than 
8 cm away from the projections of tracks into the calorimeter. Particle identification from 
dE/dx and the RICH is used on all charged particle candidates. Pions, kaons, and protons 
must be positively and uniquely identified. That is: pions must not satisfy interpretation 
as kaons or protons, and kaons and protons obey similar requirements. Charged particles 
must not be identified as electrons using criteria based on momentum, calorimeter energy 
deposition, and dE/dx. 

The invariant mass of the decay products from the following particles must lie within 
limits determined from MC studies: 7r° (M^^ = 120 - 150 MeV), rj (M^^ = 500 - 580 MeV), 
r/ (M^+^-^o = 530-565 MeV), uo (M^+^-^o = 740-820 MeV [M^+^-^o = 760-800 MeV for 
the ujpp final state]), {Mk+k~ = 1-00 - 1.04 GeV), and A (Mp^- = 1.1136 - 1.1180 GeV). 
For 7r° — i> 77 and 77 — >■ 77 candidates in events with more than two photons, the combination 
giving a mass closest to the known vr*^ or t] mass is chosen, and a kinematically constrained 
fit to the known parent mass is used. Fake vr^'s and ri's are suppressed with lateral electro- 
magnetic shower profile restrictions. For 77 tt+tt^tt'^ and u ti^ti^ti^, the vr^ is selected 
as described above, and then combined with all possible combinations of two oppositely 
charged pions choosing the combination that is closest to the 77(0;) mass. A kinematically 
constrained fit is used for neither of these modes, nor for (p's or for A's. For A — ^ pvr^, a fit 
of the pTT~ trajectories to a common vertex separated from the e~^e~ interaction ellipsoid is 
made. Contamination from Ks decays is eliminated primarily by the energy and momentum 
requirements imposed on the event, and by particle identification. 

Energy and momentum conservation requirements are imposed on the reconstructed final 
state hadrons, which have momentum pi and combined measured energy -Evis- We require the 
measured scaled energy E^is/Ecm be consistent with unity within experimental resolution, 
which varies by final state. We require |Spi|/ii^cm < 0.02. Together these requirements 
suppress backgrounds with missing energy or incorrect mass assignments. The experimental 
resolutions are smaller than 1% in scaled energy and 2% in scaled momentum. 

For the final states 2(7r+7r^)7r°, cuvr+vr", and ujK~^K~, an additional cut is applied to 
remove a background of radiative events. When the photon is combined with a low-energy 
photon signal, it can imitate a vr". We require (-E4tracks + E^)/Ecm < 0.995, where E^ is the 
energy of the highest energy photon. 

In order to compute Qh in modes with two or more charged pions, two ir^'s or an 77, 
it is necessary to remove the contribution from the transitions ipi^S) J/ipX, where 
X = TT+TT", 7r°7r°, or i] l3|. Accordingly we reject events in which the mass of any of the 



following falls within the range 3.05 < m < 3.15 GeV: the two highest momentum oppositely 
charged tracks, the recoil mass against the two lowest momentum oppositely charged tracks, 
or the mass recoiling against the 27r'''s or rj. 

For every final state, a signal selection range in E^^s/Ecm is determined by Monte Carlo 
simulation, and a sideband selection range is defined to measure background. Final states 
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with the intermediate t], uj, or particles must satisfy a scaled energy signal selection range 
requirement identical to the corresponding mode without the intermediate particle. For 
example, the scaled energy signal selection range is the same for (f)K^K~ and K^K~K^K~ . 
For final states with an 77, u, or 0, the event yield is determined from signal and sideband 
selection ranges of the intermediate particle mass. The scaled energy and mass signal and 
sideband selection ranges for ip{2S) and continuum data are identical, and are listed by 
mode in Tables HI and ITTl 

TABLE I: Signal and sideband selection ranges for the scaled energy, Ey[s/Ecm, by mode. 



mode 


signal 


sideband 


2(7r+7r-) 


0.98-1.02 


0.96- 


-0.98, 


1.02- 


1 


04 


2(7r+7r-)7r° 


0.98-1.02 


0.96- 


-0.98, 


1.02- 


1 


04 




0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 


2{K+K~) 


0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 




0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 


ppK^K~ 


0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 


AA 


0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 


AAtt+vt" 


0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 


ApK+ 


0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 


ApK^TT'^Tr~ 


0.99-1.01 


0.98- 


-0.99, 


1.01- 


1 


02 



TABLE II: The mass signal and sideband selection ranges for modes with an r], uj, or <j). Modes 
with these particles must satisfy a scaled energy signal selection requirement identical to the cor- 
responding mode without the intermediate particle. 



mode 


Variable 


signal (GeV) 


sideband (GeV) 


UJTT^TT^ 




0.74-0.82 


0.70-0.74, 0.82-0.86 


LuK+K- 




0.74-0.82 


0.70-0.74, 0.82-0.86 


LOpp 




0.76-0.80 


0.74-0.76, 0.80-0.82 






0.50-0.58 


0.455-0.50,0.58-0.615 






0.530-0.565 


0.51-0.53,0.565-0.58 


rj'Sn 


^7r7r77 


0.945-0.970 


0.9325-0.945,0.97-0.9825 




rriK+K- 


1.00-1.04 


1.04-1.08 




rriK+K- 


1.00-1.04 


1.04-1.08 


(ppp 


rriK+K- 


1.00-1.04 


1.04-1.08 
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Figures 2] and 121 show the distributions of -EvisZ-E-cm and resonance mass, where relevant, 
for all final states. The V'(2S') data is represented by the points with error bars, the open 
histogram is signal Monte Carlo with arbitrary normalization, and the shaded histogram is 
the scaled continuum. The arrows indicate the scaled energy and resonance mass selection 
ranges. The requirements of energy and momentum conservation together with particle 
identification result in small backgrounds. The dominant background in most modes at the 
ip{2S) is from the continuum. 

Event totals are shown for both the ip{2S) and the continuum in Table ITTTl where S^(^2S) 
(Sco) is the number of events in the signal region and -B^(2S) (-Bco) the number of events 
in the sideband region in ip{2S) (continuum) data. Under the assumption that interference 
between ip{2S) decay and continuum production of the same final state is absent, the number 
of events attributable to each ip{2S) decay mode, Ns, is 

Ns = 5*^(25) - -8^(25) - fiSco - -Bco), (2) 

where / is mode dependent and listed in Table 11111 We observe a statistically significant 
signal in all modes except (ppp. The signal with the least statistical significance is ujpp (i.Ja). 

The efficiency, e, for each final state is the average obtained from MC simulations for 
both detector configurations; the two values are typically within a few percent (relative) to 
each other. No initial state radiation is included in the Monte Carlo, but final state radiation 
is accounted for. The efficiencies in Table UTTl include the branching ratios for intermediate 
final states. 

We correct A^^ for the efficiency and normalize to the number of 1^(23) n^n^ J/ip, 
J/ip —>■ fi~^fi~ decays in the data, which has been determined previously to be (6.75 ± 
0.12) X 10^, where the error accounts for data and MC statistics (0.7%) and for uncertainties 
in trigger efficiency (0.5%) and angular distributions (1.5%). Complete details may be 
found in (6|. The resulting relative branching ratios are listed in Table UTTl The absolute 
branching ratios are determined using B{%l){2S) 7i^7i~J/ip) = 0.323 ± 0.013 [l| and 
-) = (5.88 ±0.10)% Q. Qh values are determined using the absolute ip{2S) 
branching ratios determined in this analysis and J/ip branching ratios from I*]. 

The systematic errors on the ratio of branching fractions share common contributions 
from the number of 4'{2S) 7i^7i~J/ip, J/ip decays (1.7%), uncertainty in / 

(1.5%), trigger efficiency (1%), electron veto (1%), and Monte Carlo statistics (2%). Other 
sources vary by channel. We include the following contributions for detector performance 
modeling quality: charged particle tracking (3% per track), vr^/r/ finding (4.4%), A finding 
(3%), Tc/K/p identification (1%/ identified tt/K/p), and scaled energy and mass resolutions 
(2%). The systematic error associated with the uncertainty in the level of background 
is obtained by recomputing the branching ratio when the background at the ip{2S) and 
the continuum are coherently increased by la. Since the background in many modes is 
small, the Poisson probability for the observed number of background events to fluctuate 
up to the 68% C.L. value is calculated and interpreted as the uncertainty in the level of 
background. The efficiencies were determined using events generated according to phase 
space. Many of the modes studied may have resonant submodes. A preliminary analysis 
finds the following prominent examples: 2(7? "•"vr") where pmr is dominant, cuvrvr where biTi 
is dominant and ci;/2(1270) is significant, K~^K~7T~^7r~ where K*{892)K7t and pK^K~ are 
large and ^vr+vr" is small, and rj'iTT where rjir^p^, rjp^TT^ , and ?7p^7r+ dominate. Allowing for 
the presence of resonant submodes changes the efficiency by less than 5% relative to the non- 
resonant efficiency for each mode we have studied. To account for this and the possibility of 
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FIG. 1: Distributions of the scaled visible energy E^is/Ecm for labeled final states. In each figure, 
entries are shown for ip{2S) data (points with error bars), signal Monte Carlo with arbitrary 
normalization (open histogram), and scaled continuum (shaded histogram). The vertical arrows 
indicate ends of signal selection ranges. 



7 




m(37T) (GeV) m(K*K^) (GeV) 01(717177) (GeV) 

FIG. 2: Distributions of the intermediate particle mass for labeled final states. In each figure, 
entries are shown for ip{2S) data (points with error bars), signal Monte Carlo with arbitrary 
normalization (open histogram) and scaled continuum (shaded histogram). The vertical arrows 
indicate ends of signal selection ranges. 

resonance polarization, we assign a 10% modeling systematic error to all modes. Systematic 
uncertainties are significant for all modes and the dominant error for many. 

Figure 01 shows the branching ratios measured in this analysis and a comparison to pre- 
vious measurements . In Table IIVI the values of Qh for those final states with measured 
branching ratios at the J/'?/' 0] are given. Figure ID shows the Qh values from Table IIVI and 
a comparison to previous measurements. This is the first measurement of Qapk- The values 
of Qh appear to be independent of the final state with no significant differences between 
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TABLE III: For each final state h: the number of events in the signal region, Sco, and background 
from sidebands, Bco, in continuum data; the scale factor, /; the number of events in the signal 
region, -S'^(25), and background from sidebands, -6^(25), in ijj{2S) data; the number of events 
attributable to ip{2S) decay, Ns, computed according to Equation (2); the average efficiency, 
e; the ratio B{^{2S) h) /B{ip{2S) — > 7r^7r~ J/i/j, J/ip and the absolute branching 

fraction with statistical (68% C.L.) and systematic errors. For t/Stt, the two decays modes 77 — > 77 
and rj ^ Sir are combined on line t/Stt. The branching ratios reported here explicitly exclude 
contributions from il){2S) — J/il^X, where X = tt+tt", tt^tt'^, or rj. 



mode 

h 


Q 

>-'co 


D 

-Deo 


J 


^tp{2S) 




AT 


e 


B{ip(2S)^h) 

7r~l~7r J / Tp 

(10-4) 


oyipylD ) — > 11) 
(10-4) 


2(7r+7r-) 


1437 


39 


0.2680 


700 


19 


306.3 


0.4317 


105.1 ± 11.5 ± 13.6 


2.0 ±0.2 ±0.3 


2(7r+7r-)7r° 


343 


28 


0.2549 


1742 


39 


1622.7 


0.1927 


1247.1 ± 46.2 ± 171.3 


23.7 ±0.6 ±3.3 




45 


7 


0.2412 


426 


38 


378.8 


0.1338 


419.5 ±33.7 ±58.5 


8.0 ±0.5 ±1.1 




28 


2 


0.2507 


271 


49 


215.5 


0.0595 


536.8 ±62.8 ±92.6 


10.2 ± 0.8 ± 1.8 




19 


10 


0.2676 


50 


1 


46.6 


0.0186 


370.2 ±66.8 ±59.1 


7.0 ± 1.1 ± 1.3 




















8.5 ± 0.7 ± 1.6 


rj'Sir 








0.2349 


14 


2 


12.0 


0.0079 


225.1 ± 113.9 ±69.9 


4.3 ± 1.5 ± 1.2 


K+K-TT+TT- 


825 


103 


0.2694 


1018 


49 


774.5 


0.3372 


340.2 ± 20.6 ± 44.0 


6.5 ±0.3 ±0.8 




26 


12 


0.2738 


72 


18 


50.2 


0.1545 


48.1 ±15.7 ±9.9 


0.9 ±0.2 ±0.2 


ujK+K- 


57 


11 


0.2091 


92 


8 


74.4 


0.1110 


99.3 ± 17.8 ± 14.8 


1.9 ±0.3 ±0.3 


2{K+K-) 


100 


10 


0.2727 


81 


2 


54.5 


0.2605 


31.0 ±8.4 ±4.6 


0.6 ±0.1 ±0.1 


(j)K+K- 


45 


16 


0.2729 


46 


4 


34.1 


0.1331 


37.9 ± 12.2 ±6.8 


0.7 ±0.2 ±0.1 




328 


31 


0.2552 


963 


31 


856.2 


0.4481 


283.0 ± 14.5 ± 36.5 


5.4 ±0.2 ±0.7 


LOpp 


18 


5 


0.2314 


29 


8 


18.0 


0.1060 


25.1 ±11.4 ±6.8 


0.5 ±0.2 ±0.1 


ppK+K- 


17 





0.2563 


29 


1 


23.6 


0.3445 


10.2 ± 3.2 ± 1.7 


0.2 ±0.1 ±0.1 


4>PP 


2 


1 


0.2472 


6 


2 


3.8 


0.1242 




< 0.18 (90%CL) 


AA 


1 





0.2431 


112 





111.8 


0.1035 


160.0 ± 19.9 ± 22.8 


3.0 ±0.3 ±0.4 


AAtt+tt" 


9 


2 


0.1775 


46 


3 


41.8 


0.0437 


141.4 ±27.1 ±26.6 


2.7 ±0.5 ±0.6 


kpK+ 


38 


5 


0.2258 


51 


8 


35.5 


0.1429 


36.8 ±9.9 ±6.3 


0.7 ±0.2 ±0.1 




17 


1 


0.1370 


26 


2 


21.8 


0.0514 


62.8 ± 23.6 ± 15.0 


1.2 ±0.3 ±0.3 



mesonic and baryonic modes. The data indicates that Qh values are in general lower than 
the current leptonic ratio. 

In summary, we have presented preliminary branching fractions for seven new decay 
modes of the iIj{2S), namely fjSir, rj'Sn, 2{K'^K'^), ppK^K^ , AAtt+tt^, ApK^, ApK^n^'n-^ , 
and more precise measurements of nine previously measured modes, which are 2(7r"'"7r-), 
2(7r"'"7r-)7r°, uj7r'^7r~, K'^K'n'^n', ^tt+tt-, ujK^K~ , (f)K^K~ , ppTT~^7r~, and AA. We also 
measure cupp and obtain an improved upper limit for (ppp. Results are compared, where 
possible, with the corresponding J/ip branching ratios to test the 12% rule. This analysis 
is part of a comprehensive study of ip{2S) multi-body hadronic decays and their resonant 
substructure and is closely related to searches for non-DD decay modes of the ^(3770). 
Further results will be presented in the near future. 
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FIG. 3: A comparison of the measurements of the ip{2S) branching ratios for the labeled final 
states from this work (circles) and, when available, from the PDG jl| (triangles). The error bars 
are the quadrature sum of the statistical and systematic errors. 
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TABLE IV: For each final state h: the branching ratio from Table III, the corresponding branching 
ratio at the J/ip from 1], and Qh from Equation (1). The ip{2S) branching ratios and Qh exphcitly 
exclude contributions from tp{2S) — > J/ipX where X = vr+vr^, vr'^vr'^, or rj. 



mode 


^(-0(25) ^ h) (10-^) 


Br(J/ih) (10^^) 


Oh (%) 


h 


this work 


PDG 


PDG 


this work 


2(7r+7r") 


2.0 ± 0.3 


4.5 ± 1.0 


40 ± 10 


5.0 ±1.5 


2(7r+7r")7r° 


23 7 ± 3 3 


30 ± 8 


337 ± 26 


7.0 ±1.1 




8.0 ± 1.2 


4.8 ± 0.9 


72 ± 10 


11.1 ±2.3 


rjSTT 


8.5 ± 1.0 










4.3 lb 2.0 










6.5 ± 0.9 


16 ± 4 


72 ± 23 


9.0 ±3.1 




0.9 ± 0.2 


1 50 ± 28 


8.0 ± 1.2 


11.4 ±3.5 




1.9 ± 0.4 


1.5 ± 0.4 


19 ± 4 


9.9 ±2.9 


2(K+K-) 


0.6 ± 0.1 








(pK+K- 


0.7 ± 0.2 


0.60 ±0.22 


15.4 ± 2.1 


4.7 ±1.6 




5.4 ± 0.7 


8.0 ± 2.0 


60 ± 5 


9.0 ±1.4 


ujpp 


0.5 ± 0.2 


0.80 ±0.32 


13 ± 2.5 


3.7 ±1.9 




0.2 ± 0.1 








(/)pp 


< 0.18(90%CL) 


< 0.26 


0.45 ± 0.15 




AA 


3.0 ± 0.5 


1.81 ±0.34 


13 ± 1.2 


23.4 ±4.6 


AAtt+tt^ 


2.7 ± 0.8 








ApK+ 


0.7 ± 0.2 




8.9 ± 1.6 


7.9 ±2.8 


Api^^vr^vr^ 


1.2 ± 0.4 
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